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Abstract We propose a quasi-Monte Carlo algorithm for pricing knock-out and 
knock-in barrier options under the Heston (1993) stochastic volatility model. This 
is done by modifying the LT method from Imai and Tan (2006) for the Heston model 
such that the first uniform variable does not influence the stochastic volatility path 
and then conditionally modifying its marginals to fulfill the barrier condition(s). 
We show this method is unbiased and never does worse than the unconditional al- 
gorithm. Additionally the conditioning is combined with a root finding method to 
also force positive payouts. The effectiveness of this method is shown by extensive 
numerical results. 



1 Introduction 

It is well known that the quasi-Monte Carlo method in combination with a good 
path construction method, like the LT method from Imai and Tan [9], can be a help- 
ful tool in option pricing, see, e.g., [4, 11]. The integrand functions usually take the 
form max(/,0) and a good path construction will somehow align the discontinu- 
ity of the derivative along the axes. However, as soon as other discontinuities, in 
the form of barrier conditions, are introduced, the performance of the quasi-Monte 
Carlo method degrades. This is also the case for the Monte Carlo method for which 
in [6] a conditional sampling method has been introduced to alleviate this problem. 

In previous work [1, Version 2] we have introduced a conditional sampling 
method to deal with barrier conditions in the Black-Scholes setting that can be used 
in combination with the LT method. In that paper we have shown that such a scheme 
always performs better than the unconditional method. Here we consider the more 
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realistic Heston model [8], which has a stochastic volatiUty component, and derive 
an algorithm to do conditional sampling on barrier conditions under this model. 



2 The LT method for Heston under log prices 

Assume a Heston world [8] in which the risk-neutral dynamics of the asset are given 
by 

dS{t) = rS{t)dt+^/v(t)S{t)dw\t), 
dV{t) = {e - V{t))Kdt + a^/v{t)dW^{t), 
dW\t)dW^{t)= pdt, 

where S{t) denotes the price of the asset at time f, r is the risk-free interest rate, k 
is the mean-reversion parameter of the volatility process, G is the long run average 
price variance and a is the volatiUty of the volatiUty. We assume the Feller condition 
2k9 > <J^ such that the process V{t) is strictly positive. The parameter p controls 
the correlation between the log-returns and the volatility. A useful observation is 
that one can write 

w\t) = pW^{t) + p^w\t), 

where W^{t) and W^{t) are independent Brownian motions. This corresponds to the 
Cholesky decomposition of the correlation structure. When resorting to Monte Carlo 
techniques for pricing options under this model, asset paths need to be discretized. 
For simpUcity we assume that time is discretized using m equidistant time steps At = 
T /m, but all results can be extended to the more general case. The notations Sj^ and 
Vk will be used for S{kAt) and V{kAt), respectively. We use the Euler-Maruyama 
scheme [10] to discretize the asset paths in log-space (see also [5, Sect. 6.5] w.r.t. 
transformations of variables) and sample the independent Brownian motions and 
W^hy using independent standard normal variables and Z^;fork = 0,...,m—l, 

log5ft+i = logSk +(r-^^At + ^/9kVA't (^pZ^+j + ^l-p^Zf^j^ , (1) 

Vk+i=Vk + ie-Vk)KAt + a^kVA'tZl^i. (2) 

For our method it is important that V is sampled solely from Z^ and to switch to 
log-space. This will be explained in the next sections. 

Write Z = (Zj ,Z^,Z2 jZj , . . . ,Z^)' G S?"\ where the prime is used to denote the 
transpose of a vector. Then Z is multivariate standard normally distributed. Assum- 
ing a European option payoff represented as 



max(/(Z),0) 
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one usually simulates the function /(Z) by mapping a uniform variate u in the unit 
cube to z by applying the inverse cumulative distribution function <P^. The method 
of then identifying Z with z will in this paper be called the standard Monte Carlo 
method (MC). When using quasi-Monte Carlo (QMC), the uniform variates are re- 
placed by a low-discrepancy point set. Our conditional sampling scheme will use the 
influence of the first uniform variable ui to try and force the barrier conditions to be 
met. For this we will employ the LT method which will use the mapping Z = Qz for 
a carefully chosen orthogonal matrix Q. This means that in (1) and (2) we take, for 
k = 0,...,m—l. 



2m 2m 
n=l 



^l+l = ^2k+l,nZn and Zf^i = l2k+2,nZn, 



where qi j denotes the element from the matrix Q at row / and column j. We remark 
that, for ease of notation, we will write /(Z), f{z),f{u) or f{Si,. . to denote 
the function / from above in terms of normal variates Z or z, uniform variates u or 
just the stock path S\,...,S,n- 

In what follows the notation Q.^- denotes the kth column of Q and g^. denotes the 
kth row. The LT method chooses the matrix Q according to the following optimiza- 
tion problem: 

maximize variance contribution of / due to Mi dimension 

Q-k 



subjectto = 

{Q*j,Q-k)=0, j = l,...,k-l, 

where Q* denotes the columns of Q that have already been optimized in the pre- 
vious iterations. The algorithm is carried out iteratively for k = 1,2, . . .,2m so that 
in the Mi optimization step the objective function ensures that, given colunms Q*j, 
j = l,...,k—l which have already been optimally determined in the previous iter- 
ations, the variance contribution due to the kth dimension is maximized while the 
constraints ensure orthogonality. Being able to express the variance contribution for 
each component analytically for general payoff functions / can be quite compli- 
cated. Therefore, Imai and Tan [9] propose to approximate the objective function by 
linearizing it using a first-order Taylor expansion for z = z 4- Az, 



2m 



df 



t=l f^Zk 



Mk- 



Using this expansion, the variance contributed due to the Mi component is 

dzk 



df 



z=t 



4 



Nico Achtsis, Ronald Cools, and Dirk Nuyens 



The expansion points are chosen as % = (1 , . . . , 1 , 0, . . . , 0), the vector with k—l 
leading ones. The optimization problem becomes 



maximize I -r — (3) 
subject to \\Q-k\\ = 1) 

{Q*j,Q-k)=o, j = i,...,k-i. 

The original Imai and Tan paper [9] considers a European call option to illustrate 
the computational advantage of the LT method under the Heston model. In their 
paper the stochastic volatility is described in [9, Sect. 4.2] and we will revisit their 
method in Section 4. For ease of illustration we also consider the payoff function 
inside the max-function to be that of a European call option 

f{z)=S,n-K 

where K is the strike price. For notational ease, we introduce the following func- 
tions: 



■^vVk 
ft = ayffkVAi. 

Note that all the above functions /' depend on Z. Similar to |9|, to find the partial 
derivatives dSm/dzi needed for the optimization algorithm, we obtain the recursive 
relations (with initial conditions dlogSo/ dzi = and dVo/dzi = 0) 

aiog5<;+i _ dlogSk , dVk 1 



dzi ^ dzi - + -d^fk + \Pl2k+i,i+^i-p^q2k+2,ijfk, (4) 
where k goes from to m — 1. The chain rule is used to obtain 

dSm ^ logoff, 

dzi "* dzi 

We will use the following lemma to calculate the optimal transformation matrix. 



Lemma 1. The recursion 



Conditional sampling for barrier option pricing under the Heston model 5 

Fk+i =akFk + hqk, 

Gk+i = CkGk + d^qk + e^Fk, 

with initial values Fq = Go = can be written at index k+ 1 as a linear combination 
oftheq£,i = 0,...,k, as follows 

k k 

Fk+1 = n 

e=o j=e+i 

k / k k k t-l \ 

Gk+l =Y,^i\'^i n ^j+^i 1^ n n • 

e=o \ j=t+i t=e+i v=t+i v=e+i J 

Proof. The formula for F^+i follows immediately by induction. For the expansion 
of Gk+l we first rewrite this formula in a more exphcit recursive form 

k k k-l k k t-l 

Gk+i=Y,ii'^( n 1^ n n 

^=0 j=t+i t=o t=t+i v=t+i v=e+i 

k k k k /t-l f-1 \ 

= Y.^i'^i n '^j+T,^' n L^^^^ n ■ 

i=o j=t+\ t=i v=t+i \e=o v=e+i J 

The part in-between the braces equals F, and the proof now follows by induction on 
k. □ 

A similar result is obtained if the second recursion is replaced by Gk+i = c^Gk + 
dkqk + d'kq'k + SkFk- Furthermore the coefficients in the expansion for qi and q'f^ can 
cheaply be calculated recursively. Using this lemma, we can make the log-LT con- 
struction for the Heston model explicit in the following lemma. 

Proposition 1. The column vector Q k that solves the optimization problem {3) for 
a call option under the Heston model is given by Q.k = ±v/|l v|| where 

m-l t-l 
V2e+l= Smfh+Sn^ft Y ft n /v' 

t=e+i v=e+i 

for i = 0,...,m—l. 

Proof. By [9, Theorem 1] the solution to the optimization problem (3) is given by 
where v is determined from 

, _ dSm _ ^ dlogSm 
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With the help of Lemma 1 we find from (4) and (5) 

^=0 V t=f.+i v=e+i ) ^=0 

from which the result now follows. □ 

Note that since Sm and all functions /' depend on Z, the vector v changes in each 
iteration step of (3) as the reference point z is changed. 

This construction can also be used for a put option with payoff 

/(z)=^-5„. 

In case of an arithmetic Asian option, the payoff is given by 

1 m 

.7=1 

In that case the optimization problem (3) contains the sum of partial derivatives 



dzk 



Z=Zk 



1 y ^ 



Z=Zk 



It is thus straightforward to use the results for the call option in Proposition 1 to 

construct the optimal transformation matrix for the arithmetic Asian option. 

Crucial to our conditional sampling algorithm is that we modify the LT construc- 
tion by forcing all odd elements in the first column of Q to zero, i.e., q2k+i,i = for 
k ~0, . . . ,m — I, while retaining the orthonormality of the columns. This effectively 
removes the influence of zi to Zj. and thus % for all k. In the next section we will 
show this leads to an elegant conditional sampling scheme. 

Corollary 1. Under the condition that q2i+i,i = Ofor i = 0,...,m—lwe have that 

the elements q2i+2,\ M have the same sign. 

Proof. From Proposition 1, for A: = 1, we find that q2i+2,i is proportional to V2f+2, 
i.e.. 



V2^+2 = Smy^e\/Ai^l-p^, 



which is always positive, and q2e+\ = V2e+i = 0. Proposition 1 now shows that the 
optimal choice is given by ±v/|| v|| from which the result follows. □ 



3 Conditional sampling on log-LT 



For expository reasons assume for now an up-&-out option with barrier B, 
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=max(/(5i,...,5m),0) l|max5A- <b| . (6) 

The condition at time that the asset stays below the barrier can then be written, 
for A: = 0,. . . ,m — 1, as 



k ^2 

= logSo + r{k+l)At- At 4- 

2m 

Zn 



k I — Zm / I 

+ £yy^v^£ ipq2e+i,n + Ji-p^q2e+2,„ 

e=0 n=2 V 

+ ZlVAi^Jl-p^Y, \f^eq2e+i,i 



< logs, 

where we have used ^2^+1,1 = 0- For notational ease we define the function 



, logB/So-r(k+l)At + Atl^l,Vf/2 

ik[o,Z2:2m) = , . —7 1=- '■ 

VAt^i - p^L}^o VVme+2,1 



iLo V^V^f I^^j (^pq2i+Ln + \/l-p2^2£+2,«) Z„ 

v^\/l-P^lLo V^eq2e+2,i 



(7) 



Here the notation Z2:2m is used to indicate the dependency on 22, . . . ,Z2m. but not zi- 
Note that depends on all other market parameters as well, but this dependency is 
supressed not to clutter the formulas. Because of the assumption that q2k+l,i = for 
all k, V can be sampled independently of zi. This means the barrier condition can 
be written as a single condition on Z]_, i.e., 



Zl < nmiIl(B,Z2:2m) if \f^eq2e+2,i > 0, 



Zl > maxIl(B,Z2:2m) if Y \/^^2^+2,l < ^• 



k 

L 

and 

k 
3) 



The condition on zi was here derived for an up-&-out option for ease of exposi- 
tion. The modifications for more complex barriers can easily be obtained from here. 
Table 1 gives an overview of the conditions on zi for the basic barrier types and 
shows that these conditions can easily be combined for more complex types. 
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Type 


£*=0 \/^?2<+2,l > 


U&O (B) 
D&O (B) 

U&l (B) 

D&l (B) 
U&O + D&O (B, > B2) 
U&O + D&l (Bi > B2) 


Zl < mmtrk{B,Z2:2m) 
Zl >maX^-Il{B,Z2:2m) 
Zl > min^I](B,Z2:2m) 
Zl < maXj-Il(B,Z2:2m) 
Zl e {maXti;f(B2,Z2:2m),min*r4(Bi,Z2:2„,)) 

Zl < mm{max*-r^(B2,Z2:2m),min*-Jl{Bi,Z2:2m)} 


Type 


ELo \/W'?2/'+2.1 <0 


U&O (B) 

D&O (B) 

U&l (B) 

D&l (B) 
U&O + D&O (Bi > B2) 
U&O + D&l (Bi > B2) 


Zl > maX^-Il(B,Z2:2m) 
Zl < mintIl{B,Z2:2m) 
Zl <maXtIl(B,Z2:2m) 
Zl >niin;fcI]t(B,Z2:2m) 
Zl e (maXiIi(Bi,Z2:2m),ininiJi(B2,Z2:2m)) 
Zl > max{maXi.rt(Bl,Z2:2m),min;(:ri(B2,Z2:2m)} 



Table 1 The barrier constraints on zi for different types of barriers: up-&-out (U&O), down-&-out 
(D&O), up-&-in (U&l), down-&-in (D&l) and some combinations. 



We now show the main results on our conditional sampling scheme. Again, for 
expository reasons, specialized for the case of the up-&-out option from above. This 
result can easily be modified for other payout structures in the same spirit as the 
results in Table 1 . The following theorem holds for both the Monte Carlo method as 
for a randomly shifted quasi-Monte Carlo rule. 

Theorem 1. For the up-&-out option (6) and assuming that we fixed qu+i^i > for 
i = 0,...,m—l {see Corollary 1) the approximation based on sampling 

8{zi,.--,Zm) = ^ ^imnIX(£,Z2:2m)^ max (/(zi ,Z2, • • • ,Zm) , 0) 

where, using the relation zi = <P~^{ui), 

Zl = <P~^ (ui nunIX(5,Z2:2m)^ , (8) 
is unbiased. Furthermore, if we denote the respective unconditional method by 
g{z\,...,Zm)=rQ3ix{f{Si,...,Sm),0) l|max4<5|, 

where the S\, Sm are obtained directly from z\,. . . ,Zm without using (8), then, 
when using the Monte Carlo method or a randomly shifted quasi-Monte Carlo 
method, the conditional sampling has reduced variance, i.e., Var[^] < Var[g]. Fur- 
thermore the inequality is strict i/Pfmaxi^jt >B\>Q and E[g] > 0, i.e., if there is 
any chance of knock-out and positive payoff. 



Proof. The proof can be constructed similar to [1, Theorem 3, 4 & 5] from our 
previous work. □ 
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The previous result shows that the proposed conditional algorithm can never do 
worse than its unconditional variant. Furthermore, the more chance there is on a 
knock-out the more effect the conditional algorithm will have. This can be observed 
in the examples in Section 5. 

Remark. The conditional sampUng was appUed to zi (or, equivalently, to mi) to keep 
the asset from knocking out (or in). Taking it one step further one could try to add 
an additional bound on zi, keeping Zi-.im constant, in order to force a strictly pos- 
itive payout. This is more involved than the barrier condition however as for more 
complicated payoffs than calls and puts there might not exist analytical formulae 
such as in Table 1 to condition zi . It is interesting to note that for calls and puts the 
same formulas can be used as in Table 1, only now restricting the functions to 
rm{K,Z2:2m)- Adding this constraint to the existing barrier conditions is straightfor- 
ward. Root finding methods can be employed for more complex payout structures. 
See our previous work [1] for a detailed analysis of root finding for Asian options. 



4 The original LT method for Heston 

We mentioned previously that it is essential for our method to switch to log prices. 
To illustrate the problem, we introduce the LT method for the Heston model as in [9] 
in this section. Also here we can obtain an explicit form of the orthogonal matrix Q 
in the spirit of Proposition 1, but the conditional sampUng scheme from the previous 
section is not applicable as we will show. The Euler-Maruyama discretizations for 
S{t) and V{t) are given by 



compare with (1) and (2). For ease of notation, we introduce the following func- 
tions: 





1 - KAt + 



■Z, 





10 Nico Achtsis, Ronald Cools, and Dirk Nuyens 

Note that all the above functions /' depend on Z. The recursion relations for the 
partial derivatives become 







dZi 


dzi- 






dzi 





if^^r, r,f3 , ^ /i _„2f3 



fk +^2k+\,ifh 

for A: = 0, . . . ,OT — 1, and initial conditions dSo/dzi = and dVo/dzi = 0. With this 
notation we obtain the LT construction for the Heston model in explicit form. 

Proposition 2. The column vector v = Q.^ that maximizes the optimization prob- 
lem {3) for a call option under the Heston model is given by Q.^ = ±v/||v|| where 

m— 1 m— 1 m— 1 r— 1 

V2m=/ip n fj+fi I ft n /v n 

j=e-^ 1 t=£-\- 1 v=f + 1 v=^+ 1 

m— 1 



V2i+2=f!Jl-p^ n //' 

fori = 0, ... ,m — 1. 

Proof. The proof is similar to Proposition 1 , again making use of Lemma 1 . □ 

To show the advantage for conditional sampling of the log-LT method (as ex- 
plained in Sections 2 and 3) over this version we consider again the up-&-out option 
with payoff 

g{Su---,Sm) = max(/(5i,...,5m),0) ijmax^A; <b| . 
The barrier condition at an arbitrary time step ti^+i takes the following form: 



Sk+i = Sk i^l+rAt+^JVkVAi (^pZ^+i + ^J 1 -p^Zf+i 

= •^0 0^1+ rAt + \j9f/Ai £ (^pq2i+i,n + \ - p'^qn+2,^ Z« j 

Trying to condition on zi as we did in the log-LT model leads to the following 
condition: 

k / I— , \ 5 



=0 

where 



n {Ai + \f9f/AiJ\-p'^q2i+2,\Z'^ < 



So 
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Ae = I + rAt + yVeVAt [Pq2i+l,n + y i - P^q2e+2,n) Zn, 
n=2 V / 

and we used again q2i+\.\ = to sample V independent of z\. To satisfy the condi- 
tion on zi, a fc + 1-th order polynomial must be solved in order to find the regions 
where the above condition holds. To find the global condition, one has to solve 
polynomials of degrees 1 to m, and then find the overlapping regions where all con- 
ditions hold. This quickly becomes impractical and we therefore use the log-LT 
method which does not have this drawback. 



5 Examples 

Up-&-out call and put 

Consider the up-&-out call and put options with payoffs 

Pc{,^\i • • • ; Sm) 
Pp(.^\i ■ ■ ■ i^m) 

The fixed model parameters are r = 0% and k: = 1. Furthermore, time is dis- 
cretized using m = 250 steps and thus our stochastic dimension is 250. The re- 
sults for this example are calculated using a lattice sequence (with generating vector 
exod8_base2_ml3 from [12] constructed using the algorithm in [2]). The improve- 
ments of the standard deviations w.r.t. the Monte Carlo method for different choices 
of p. So, Vq = = G, K and B are shown in Table 2. The results for the call and 
put option seem to be consistent over all choices of parameters: the new conditional 
scheme (denoted by QMC+LT+CS) improves significantly on the unconditional 
LT method (denoted by QMC+LT). Note that the QMC+LT method uses the con- 
struction of Proposition 2. Adding root finding (denoted by QMC+LT+CS+RF), 
to force a positive payout, further dramatically improves the results. The improve- 
ment of the QMC+LT+CS method for the put option is even larger than that for the 
call option. This difference should not come as a surprise: when using conditional 
sampling on a knock-out option, zi is modified such that the asset does not hit the 
barrier. In case of an up-&-out call option, the asset paths are essentially pushed 
down in order to achieve this. The payout of the call option however is an increas- 
ing function of §„, so that pushing the asset paths down has the side-effect of also 
pushing a lot of paths out of the money. For the put option the reverse is true: the 
payout is a decreasing function of §„, meaning that pushing the paths down will re- 
sult in more paths ending up in the money. Root finding can be used to control this 
off-setting effect in case of the call option, this effect is clearly visible in Table 2. 



; max 



(5„-/s:,0)l|max4<£|, 



: max {K - Sm, O) I \ maxSk < B 
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{Vo = 6 = a,p,So,K,B 



(0.2,- 
(0.2,- 
(0.2,- 
(0.2,0, 
(0.2,0, 
(0.2,0, 



0.5,90,80,100) 
0.5, 100, 100, 120) 
0.5,110,100,150) 

5,90,80,100) 
5,100, 100,125) 
5,110,100,140) 



(0.3,- 
(0.3,- 
(0.3,- 
(0.3,0, 
(0.3,0, 
(0.3,0, 



0.5,90, 80,100) 
0.5,100, 100,120) 
0.5,110,100,150) 
5,90,80,100) 
5, 100, 100, 125) 
5,110,100,140) 




Put 



0.05 
0.06 
0.77 
0.05 
0.10 
0.34 



(0.2,-0.5,90,80, 100) 


367% 


331% 


184% 


9.02 


(0.2,-0.5,100,100,105) 


279% 


235% 


126% 


7.76 


(0.2,-0.5,110,100,112) 


I 298% ^ 


263% 


123% 


4.44 


(0.2,0.5,90,80, 100) 


L 361% M 


376% J 


148% 


6.05 


(0.2,0.5,100,100,105) 




298gH 


131% 


5.33 


(0.2,0.5,110,100,112) 




325^^* 


149% 


2.98 


(0.3,-0.5,90,80,100) 


^^^383% 


348% 


137% 


10.3 


(0.3,-0.5,100,100,105) 


HP 260% 


243% 


144% 


8.65 


(0.3,-0.5,110,100,112) 


214% 


187% 


129% 


5.38 


(0.3,0.5,90,80,100) 


380% 


294% 


160% 


6.44 


(0.3,0.5,100,100,105) 


304% 


272% 


174% 


5.57 


(0.3,0.5,110,100.112) 


305% 


279% 


124% 


3.33 



Table 2 Up-&-out call and put. The reported numbers are the standard deviations of the MC 
method divided by those of the QMC+LT+CS+RF, QMC+LT+CS and QMC+LT methods. The 
MC method uses 30720 samples, while the QMC methods use 1024 samples and 30 independent 
shifts. The rightmost coltmin denotes the option value. 



These numerical results are illustrated in terms of in Figure 1 for two parameter 
choices for the call option. 



Up-&-in call 

Consider an up-&-in call option with payoff 

P(5i, . ..,§,„)= max {§„ -K,0)l jmaxS^ > b| . 

The fixed model parameters are r = 2%, k = I and a = 0.2. Again, m = 250. Here 
we use the Sobol' sequence with parameters from [7] and digital shifting [3]. The 
standard deviations for different choices of p, Sq, Vq = d, K and B are shown in 
Table 3. The improvements of the conditional scheme are extremely high for this 
case. Note the impact of the correlation on the results: the improvement for p = 0.5 
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(0.2, -0.5, no, 100, 150) (0.3, 0.5, 100, 100, 120) 




29 2'" 2" 2'2 2' 2'" 2" 2'^ 

iV iV 



Fig. 1 Up-&-out call convergence plots for two options with different parameters. The fixed pa- 
rameters are r = 0% and k = 1. The different choices for (Vb = = a,p,So,K,B) are denoted 
above the figures. 



{Vo = e,p,So,K,B) 


QMC+LT+CS+RF 


QMC+LT+CS 


QMC+LT 


Value 


(0.1,-0.5,90,80,160) 


2158% 


1515% 


242% 


5.47 


(0.1,-0.5,100,100,180) 


■^2377% ^ 


1542% 


240% 


5.05 


(0.1,-0.5,110,120,200) 


^^572% M 


1545% J 


250% 


4.74 


(0.1,0.5,90,80,160) 
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Table 3 Up-&-in call. The reported numbers are the standard deviations of the MC method di- 
vided by those of the QMC+LT+CS+RF, QMC+LT+CS and QMC+LT methods. The MC method 
uses 30720 samples, while the QMC+LT+CS+RF, QMC+LT+CS and QMC+LT methods use 1024 
samples and 30 independent shifts. The rightmost column denotes the option value. 



is even approximately twice that for p = —0.5. All parameter choices indicate that 

conditional sampling on the barrier condition greatly improves accuracy. Adding 
the additional condition of the payout itself (root finding) provides another serious 
reduction in the standard deviation. 



Up-&-out Asian 

Consider an up-&-out Asian option with payoff 
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Table 4 Up-&-out Asian call. The reported numbers are the standard deviations of the MC method 
divided by those of the QMC+LT+CS+RF, QMC+LT+CS and QMC+LT methods. The MC method 
uses 30720 samples, while the QMC+LT+CS+RF, QMC+LT+CS and QMC+LT methods use 1024 
samples and 30 independent shifts. The rightmost column denotes the option value. 



(0. 1 , -0.5, 110, 120, 160) (0. 15, 0.5, 90, 80, 120) 




Fig. 2 Up-&-out Asian call convergence plots for two options with different parameters. The fixed 
parameters are r = 5%, K" = 1 and <t = 0.2. The different choices for (Vb = 0,p,So,K,B) are 
denoted above the figures. 



P{Si,. . . = max (^^ £ -K,O^I |max5* < B 

The fixed model parameters are r = 5%, K" = 1 and a = 0.2. The number of time 
steps is fixed at m = 250. We use the Sobol' sequence as in the previous example 
and the results are shown in Table 4. The results are once more very satisfactory 
with similar results as for the up-&-out call and put options in Table 2. Figure 2 
shows the convergence behaviour for two sets of parameter choices. As before, a 
significant variance reduction can be seen for our conditional sampling scheme and 
the root finding method further improves this result. 



Conditional sampling for barrier option pricing under the Heston model 

6 Conclusion and outlook 
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The conditional sampling scheme for the LT method introduced in [1] for the Black- 
Scholes model has been extended to the Heston model. This was done by consider- 
ing log prices and making the sampling of the volatility process independent of zi- 
We also obtained explicit constructions for the optimal Q matrix for the LT method 
in this case. The numerical results show that the method is very effective in reduc- 
ing variance and outperforms the LT method by a huge margin. We only considered 
an Euler-Maruyama discretization scheme for the asset and volatility processes. It 
might be interesting to see if the theory and results carry over when other simulation 
methods are used, see [13] for an overview of other methods. 
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